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Abstract 

In this paper we study the twist disclination within the elastoplastic defect theory. 
Using the stress function method, we found exact analytical solutions for all characteristic 
fields of a straight twist disclination in an infinitely extended linear isotropic medium. The 
elastic stress, elastic strain and displacement have no singularities at the disclination line. 
We found modified stress functions for the twist disclination. In addition, we calculate 
the disclination density, effective Frank vector, disclination torsion and effective Burgers 
vector of a straight twist disclination. By means of gauge theory of defects we decompose 
the elastic distortion into the translational and rotational gauge fields of the straight twist 
disclination. 

Keywords: disclinations; dislocations; gauge theory of defects; stress functions 
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1 Introduction 

Disclinations are very important and interesting lattice defects. They may be defects in warped 
and twisted materials. Disclinations have been investigated in the context of applications to 
liquid crystals as twisting discontinuities [1,2], Abrikosov lattices formed by magnetic flux 
lines in the mixed state of type-II superconductors [3] , polymers by chain kinking and twisting 
of molecules [4], Bloch wall lattices [5,6], biological structures [7], amorphous bodies [8] and 
rotation plastic deformations [9, 10]. Because disclinations cause strong elastic distortions and 
lattice bending it seems that very strong distortions are necessary in order to realize disclinations 
in crystals. 

A disclination is characterized by a closure failure of the rotation for a closed circuit round 
the disclination line. There are wedge and twist disclinations. If the Frank angle (rotation 
failure) of the disclination is a symmetry angle of lattice, then the disclination is called a perfect 
disclination. Such disclinations have been introduced by Anthony [11] and deWit [12,13]. In 
the case of a twist disclination the rotation axis is perpendicular to the disclination line. The 
smallest value of the Frank vector is 7r/2 in a cubic lattice and tt/3 in a hexagonal lattice. If 
the Frank angle is not a symmetry angle of lattice, the disclination is called partial disclination. 

'Present address: Laboratoire de Modelisation en Mecanique, Universite Pierre et Marie Curie, Tour 66, 
4 Place Jussieu, Case 162, F-75252 Paris Cedex 05, France. 
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They play an important role, e.g., in building of twin boundaries (see, e.g., [14]). Disclinations 
correspond, in general, to Volterra's distortions of the second kind (see also [15]). Thus, these 
defects are of rotational type. They are different from the so-called Prank's (spin) disclinations 
which are elementary defects in liquid crystals (see [1]). 

The traditional description of elastic fields produced by defects (e.g. dislocations, disclina- 
tions and cracks) is based on the classical theory of linear elasticity. However, classical elasticity 
breaks down near the defect line and leads to singularities. This is unfortunate since the de- 
fect core is a very important region in the theory of defects. Of course, such singularities are 
unpliysical and an improved model of defects should eliminate them. 

On the other hand, there arc other non-standard continuum models of defects, e.g., the 
nonlocal continuum model [16-21], the strain gradient elasticity [22-28] and the field theory of 
elastoplasticity which has been developed from the gauge theory of defects [29-34] . All these 
theories are successfully applied to the description of screw and edge dislocations. In this context 
the stresses have no singularities at the dislocation line. In addition, the dislocation core arises 
naturally. In particular, the field theory of elastoplasticity is a gauge theory of defects in which 
the defects cause plasticity. The corresponding gauge fields may be identified with the plastic 
distortion. By the help of this theory the elastic and plastic part of the total distortion can be 
calculated. The total distortion is defined in terms of a displacement and consists of the elastic 
and plastic part. In the case of dislocations (see, e.g., [33]) the elastic distortion is continuous 
even in the dislocation core and the plastic part becomes discontinuous. But in the case of 
disclinations the situation is less complete worked out. The stresses of straight wedge and twist 
disclinations have been calculated by Povstenko [20] in the framework of Eringcn's nonlocal 
elasticity and by Gutkin and Aifantis [26-28] by the help of strain gradient elasticity. However, 
no rotation and displacement vectors, no bend-twist and no disclination and dislocation density 
tensors were obtained in their works. In a recent paper [34] the wedge disclination has been 
investigated in the field theory of elastoplasticity. It was possible to calculate all characteristic 
field quantities. It has been seen that the disclination core may be defined quite natural in this 
framework. 

In this paper we want to extend our study for a straight twist disclination. We use the field 
theory of elastoplasticity to find nonsingular solutions for the stress and strain fields and the 
rotation and displacement fields. In addition, we investigate the relation to gauge theory of 
defects. We use the stress function method and hope to close the gap between the non-local 
and strain gradient results for the case of a straight twist disclination. In this framework we 
want to work out all geometric quantities of a twist disclination. 

2 Basic equations 

In this section we apply the field theory of elastoplasticity to the case of a straight twist 
disclination. In elastoplasticity the elastic distortion is given by [29-32] 

Pij=djUi+kj- (2.1) 

It is an additive decomposition of the elastic distortion into compatible and purely incompatible 
distortion. This decomposition can be justified by the help of the gauge theory of defects [29] . 
The displacement field Ui gives rise to a compatible distortion and the tensor is the proper 
incompatible part of the elastic distortion. 

The Burgers vector hi is defined by the help of the distortion tensor 

bi= (f (3ijdxj, (2.2) 
where 7 denotes the Burgers circuit. In elastoplasticity the linear elastic strain tensor is given 
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by means of the incompatible distortion tensor (2.1) according to 

Eij = P(ij) = ^ {diUj + djUi + Pij + fiji) , Eij = Ej,. (2.3) 

The force stress is the response quantity to elastic strain and is given by the (generalized) 
Hooke's law for an isotropic medium 

aij = 2/i ^Eij + ^_2j^ ^ijEkkJ , cTij = crji, (2.4) 

where /x, v are shear modulus and Poisson's ration, respectively. The force stress satisfies the 
force equilibrium condition 

dja^j = 0. (2.5) 

The inverse of Hooke's law reads 

Eij = ^ {^ij - Sijf^kk^ ■ (2.6) 

In the conventional disclination theory [11-13,35] the torsion tensor (linear version of Car- 
tan's torsion) is defined by 

aij ■■= e.jki{dkPii + eum^'^k) = <^jki{dk(3ii + Ciimfrnk)- (2-7) 

Anthony called it the disclination torsion (see [11]). On the other hand, it is sometimes called 
dislocation density in the theory of disclinations (see, e.g., [12, 13,35]). The y*^- was introduced 
by Mura [35] as "plastic rotation" and dcWit [12,13,36] called this quantity "disclination 
loop density". The field Pij may be identified with de Wit's "dislocation loop density". For 
a dislocation it yields ip*j = and then (2.7) has the shape of a proper dislocation density. 
Using the elastic bend- twist tensor (see, e.g., [37]) 

kij = djiJi- ifij, (2.8) 

with the rotation vector 



T^^ijkPjk, (2.9) 



Eq. (2.7) can be rewritten according to (see also [11,36,37]) 

aij = Cjki {dkEii + Cimikmk) = ^jkidkEii + Sijku - kji. (2-10) 

The index i indicates the direction of the Burgers vector, j the dislocation line direction. Thus, 
the diagonal components of represent screw dislocations, the off-diagonal components edge 
dislocations. 

The so-called disclination density tensor of a discrete disclination is defined by [11-13, 35-37] 

^ij • — ^jmn^mkin — ^jmn^m^in' (2-11) 

The index i indicates the direction of the Frank vector, j the disclination line direction. Thus, 
the diagonal components of 8^ represent wedge disclinations, the off-diagonal components twist 
disclinations. The Prank vector Oj is defined by the help of the elastic bend-twist tensor 

fl,= <J) kijdxj. (2.12) 
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Consequently, the dislocation density and the disclination density satisiy the following com- 
patibility equations (1st and 2nd Bianchi identities) 

djUij - eikiQki = 0, (2-13) 
djO.j = 0. (2.14) 

On the other hand, the theory of defects (dislocations and disclinations) can be considered 
as a gauge model of defects in solids [38,39]. The gauge group is the group ISO{3) — T(3)X) 
50(3) (r(3) - three-dimensional translational group, 5*0(3) ~ three-dimensional rotational 
group and Z) denotes the semi-direct product). In this framework, we are able to decompose 
the incompatible distortion (2.1). Namely, the incompatible distortion takes the (linearized) 
form [38-41] 

0ij = (pij + eikiWkjXi, (2.15) 

where (pij and Wij are the translational and rotational gauge fields, respectively. More precisely, 
(pij is the translational part of the generalized affine connection [42, 43] and Wij the rotational 
connection (see also [44]). The torsion and the disclination density tensor are defined by 

aij = <^jkidk4>ii + etki&kjXi, (2.16) 

^ij ^jmn^m^^in- (2*1'^) 

The disclination density tensor (2.17) is the linearized Riemann-Cartan curvature tensor or 
equivalent the corresponding Einstein tensor. It can be seen that a non-vanishing disclination 
tensor (2.17) gives a contribution to the torsion tensor (2.16). This piece may be called the 
disclination torsion. Therefore, the torsion (2.16) has a contribution from both the translational 
sector (=dislocations) and the rotational sector (=disclinations). Of course, in the case of 
dislocations (teleparallelism) , 0^ = 0, the disclination torsion is zero and only the first piece 
in (2.16), which is the proper dislocation density tensor, gives a non- vanishing contribution. If 
we compare Eq. (2.11) with (2.17), we may identify (see also [40]) 

Wij = -^\j. (2.18) 

Using Eqs. (2.15), (2.17) and (2.18), one is able to prove the equivalence between (2.10) and 
(2.16). 

The basic equation for the force stress in an isotropic medium is the following inhomogeneous 
Helmholtz equation [32] 



fl - K-'^ls\aij = aij, = ^, (2.19) 



where aij is the stress tensor obtained for the same traction boundary-value problem within the 
theory of classical elasticity. It is important to note that (2.19) agrees with the field equation for 
the stress field in Eringen's nonlocal elasticity [16, 17] and in gradient elasticity [24]. The factor 
has the physical dimension of a length and it defines, therefore, an internal characteristic 
length. If we consider the two-dimensional problem and using Green's function of the two- 
dimensional Helmholtz equation, we may solve the field equation for every component of the 
stress field (2.19) by the help of the convolution integral: 

(^iii'r)= I a{r-r')aij{r')Av{r'), (2.20) 
Jv 

with the two-dimensional Green's function 

a(r-r') = ^i^o(«(r-r')), (2.21) 
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with r = a/x^ + Here K„ is the modified Bessel function of the second kind and n = 0, 1, . . . 
denotes the order of this function. Thus, 

- K-^A)a(r) = 6{r), (2.22) 

where 6{r) := S{x)S{y) denotes the two-dimensional Dirac delta function. In this way, we 
deduce Eringen's so-called nonlocal constitutive relation for a linear homogeneous, isotropic 
solid with Green's function (2.21) as nonlocal kernel. This kernel (2.21) has its maximum at 
r = r' and describes the nonlocal interaction. Its two-dimensional volume-integral yields 



X 



a{r - r') Av{r) = 1, (2.23) 



and is the normalization condition of the nonlocal kernel. In the classical limit (k ^ ^ 0), it 
becomes the Dirac delta function 

lim a{r - r') = S{r - r'). (2.24) 

Note that Eringen [16-19] found the two-dimensional kernel (2.21) by giving the best match with 
the Born-Karman model of the atomic lattice dynamics and the atomistic dispersion curves. 
He used the choice cq = 0.39 for the length 

K;-i=eoa, (2.25) 

where a is an internal length (e.g. atomic lattice parameter) and eo is a material constant. 

Using the inverse of the generalized Hooke's law (2.6) and (2.19), we obtain an inhomoge- 
neous Helmholtz equation for every component of the strain tensor (see [32]) 

(l-K-^AjEij =Eij, (2.26) 

o 

where Eij is the classical strain tensor. Equation (2.26) is similar to the equation for the strain 
in gradient theory used by Gutkin and Aifantis [22-24] if we identify with the gradient 
coefficient (see, e.g., equation (4) in [24]). Since the strain tensor fulfils an inhomogeneous 
Helmholtz equation, we may rewrite (2.26) as a nonlocal relation for the strain 

Eij{r) = [ a{r- r') Eij{r') dv{r'), (2.27) 
Jv 

which is similar to the nonlocal relation for the stress (2.19). Thus, field theory of elastoplasticity 
may be considered as a nonlocal theory for the stress as well as the strain tensor. In contrast 
to Eringen's nonlocal theory where only the stress tensor has a nonlocal form. We assume that 
the stress and strain fields at infinity should have the same form for both the classical and 
elastoplastic field theory. 



3 Classical solution 

In this section wc present the "classical" stress field for a straight twist disclination in an 
infinitely extended isotropic body by the help of the stress function method. We assume the 
disclination line is along the z-axis and the Frank vector has the following form ft = (0, fl, 0). 
In contrast to the case of a wedge disclination or screw and edge dislocations, the situation is 
not really a two-dimensional problem for the twist disclination. In the case of a straight twist 
disclination the three-dimensional space may be considered as a product of the two-dimensional 
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xy-plane and the independent one- dimensional ^;-line [13]. In this situation the ^;-axis plays a 
peculiar role. 

The classical solution for the elastic stress fields was originally given by deWit [13] 



27r(l - I/) 



27r(l - u) r4 



'xy 



O'zv = — 



27r(l - u) 






zy 


7r(l - v) 




liQ 
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27r(l - i/) 


2 ' 




/ 



27r(l - v) 



(1 - 2!y)lnr + 



(3.1) 
(3.2) 

(3.3) 
(3.4) 
(3.5) 
(3.6) 



Obviously, the expressions (3.1)~(3.4) contain the classical singularity ~ and a logarithmic 
singularity ^ Inr in (3.6). Thus, the classical elastic stress is infinite at the disclination line. 
The reason is that the classical theory of elasticity breaks down in the disclination core so that 
in the defect core region classical elasticity fails to apply. Usually, the radius of this region is 
estimated by means of atomic models. Due to the unphysical singularities it is erroneous to 
argue that the stress has a maximum/minimum value at the defect line. 

o 

For the situation of the strain condition, E^z = 0, Eqs. (3.1)-(3.6) can be calculated by 
using the so-called stress function method in the following form 



^ d^yl -dlj -dyF ^ 

-dlj dlj dj + dz°g 

\ -dyF dj + dz°g V ) 



(3.7) 



The stress is given in terms of the stress functions /, F, g and p. In order to satisfy the force 
equilibrium the stress (Jzz has to fulfil the condition 



p = uAf = -dyg, 



(3.8) 



where A = d^^ + dyy denotes the two-dimensional Laplacian. The "classical" stress functions 
for the stress fields (3.7) are 



/ = - 

O 

F = 



lift 



9 = 



27r(l - v) 

'27r(l-i/) 
IJbQ,v 



zy In r, 
xlnr. 



7r(l - v) 



zhi r. 



(3.9) 
(3.10) 
(3.11) 
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They satisfy the following two-dimensional differential equations 

2/iOz 
(1^) 



AA/ = -,;^a,<5(r), (3.12) 



AAF = -^^^dJ{r), (3.13) 

^5=(rr^'^W- (3.14) 

o o o o 

Thus, / and F are biharmonic stress functions and g is a harmonic one. We see that F is an 

Airy stress function, / is an Airy stress function multiplied by z and on the other hand g is a 
Prandtl stress function multiplied by z (up to constant pre- factors) . 

For convenience we give the classical elastic strain of the straight twist disclination (see [13]) 



47r(l-i/) r2' 

47r(l — u) '' ' r 

which contains the "classical" singularities at r = 



(3.18) 



47r(l 

Q, xy 



= -I^7T3T:T (1 - 2-) + 72 h (3-19) 



4 Nonsingular solution 

In this section we want to consider the twist disclination in the elastoplastic field theory to 
find modified solutions without the "classical" singularities. The modified solutions are used to 
estimate the extent of disclination core, thus providing information which cannot be obtained 
by using classical elasticity theory. 

We make for the modified stress field an ansatz in terms of unknown stress functions which 
has the same form as the classical stress field (3.7) 

/ dy ~dlj -dyF \ 

'^ii=\ -dlyf dlJ d,F + d,g , (4.1) 
V -dyF d^F + d,g p J 



with the relation 



p = uAf = -dyg. (4.2) 



Substituting (4.1) and (3.7) into (2.19) we obtain three inhomogeneous Helmholtz equations 
for the unknown stress functions 
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The inhomogeneous parts of (4.3)~(4.5) are the classical stress functions. Using the same proce- 
dure as in the case of a straight edge dislocation (see [32]) in order to solve the inhomogeneous 
Helmholtz equations, we can find the solutions of (4.3)-(4.5). The solutions for the modified 
stress functions of a straight twist disclination are given by 

/ = - ^ , zy(lnr+^(l-^rKi(Atr))l, (4.6) 



27r(l - v) 

F = - if^ x{ \nr + ^[l ^ KrK^{Kr)) \, (4.7) 
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27r(l - v) 



7r(l - ly) 



— z^lnr + Ko{Kr)Y (4.8 



where the first pieces are the classical stress functions (3.9)-(3.11). 

By means of Eq. (4.1) and the stress functions (4.6)-(4.8), we are able to calculate the 
modified stress of a straight twist disclination. So we find for the elastic stress in Cartesian 
coordinates 



"2^ 



3x2) ^ __(y2 _ 3^2^ _ 2y^KrKi{Kr) - 2(y2 _ Sx^)K2{Kr) 

Ki T 

(4.9) 
(4.10) 



{(^' - y^) -^i^^- - 2y'KrJri(«r) + 2(ar2 - 3y2)jr2(Kr)|, 

(4.11) 

■{l-;^(2-«Vi^2(«r))|, (4.13) 

{l-2iy){liir + Ko{Kr)) + '^-^^^-^^^(2-K^r^K2{Kr)^ [. (4.14) 



CTzx = 



CTzv = — 



27r(l-zy) 
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27r(l - 






f 



27r(l - v) 



These stresses are plotted in Fig. 1. If we identify k = 1/^/c (c is the gradient coefficient 
used by Gutkin and Aifantis), the components of the stress (4.9)-(4.14) are in agreement with 
the stress field obtained by Gutkin and Aifantis [27, 28] in the framework of strain gradient 
elasticity by using the Fourier transform method. It is interesting to note that the stresses (4.9)- 
(4.12) caused by the straight twist disclination with the Frank vector ft = (0,f2,0) coincide 
with the stresses due to the straight edge dislocation with the Burgers vector 6 = {b, 0, 0) 
replacing flz by b (compare with equations (3.15)-(3.18) in [32]). The trace of the stress tensor 
ffefe = (Txx + (^yy + <^zz produccd by the twist disclination in an isotropic medium is 

a,, = -t^^^^{l-.rK,{.r)]. (4.15) 
7r(l — i^) L J 

We may now discuss some details of the stresses near the disclination core region in the xy-plane. 
The stresses (4.9)-(4.12) vanish at the disclination line instead of being singular as predicted 
by classical elasticity. Every component of (4.9)-(4.12) have a maximum and a minimum near 
the disclination line. Because the extreme values are of opposite sign a zero point must be at 
the defect line. In addition, the stress (4.14) has a maximum value at the disclination line. 
The extreme values may serve as a measure of critical stress level at which fracture or failure 
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can occur. Contrary to classical elasticity, stresses (4.9)-(4.14) are finite at the defect line. 
Therefore, the stress fields have no artificial singularities at the core and the maximum stress 
occurs at a short distance away from the disclination line (see Fig. 1). In fact, when r ^ 0, 
we have 



Kr 



2 



1 _ , . 1 



A'o(Kr) ^ - 7 + ln— , Ki{Kr) ^ — , K2{Kr) 



Kr' 2 (Kry 



and thus cTy 0. Here 7 denotes the Eulcr constant. It can be seen that the stresses have 
the following extreme values in the a;y-planc: |(Tj:2;(0, y)| — ^■^^^'^ 27t^{i-u) l^^l — 0.996k~^, 
\ayy{0,y)\ :^ 0.260/^5^^ at \y\ ~ 1AMk-\ Ky{x,0)\ ~ 0.260^2;^ at ~ 1.494k-i, 
|a,,(0,y)| ~ 0.399K^t^ at \y\ ^ 1.114k-i and \akk{0.y)\ ^ 0.399k ^g^^+g" at \y\ :^ 1.1Uk-\ 
The stresses axx, cryy and modified near the disclination core (0 < r < 12«; ^). The 

stress dzz and the trace akk are modified in the region: < r < 6k~^. Far from the disclination 
line (r ^ 12k^-'^) the modified and the classical solutions of the stress of a twist disclination 
coincide. In addition, it can be seen that at 2 = the stresses (4.9)-(4.12) are zero. The stress 
(Tzy has at r = the maximum value: <Tzy{0) ~ 2ir(i^iy) [(^ ^ 2;^) (7 + In -j) — i] and with = 0.3: 
azyix,0) ~ 2Tr(i-u) [Q-^ ^ ~ 0.546] (see Fig. le where a constant term proportional to Inn is 
dropped out). Consequently, one can equate the maximum shear stresses to the cohesive shear 
stresses to obtain conditions to produce a disclination of single atomic distance. 

Due to the two-dimensional symmetry it is convenient to express the stresses in cylindrical 
coordinates. The stress tensor has the following form in cylindrical coordinates 

{l-^ + 2i^2(/.r)|, (4.16) 

l-;^+2^2(Kr)|, (4.17) 

1^1 + ^ - 2K2{Kr) - 2KrKi{Kr)Y (4.18) 

l-KrKi{Kr)Y (4.19) 



CTzz = — 
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27r(l- 


v) r 
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27r(l - V 


) r I 




zsini^ 


27r(l- 


v) r 




zsintp f 



7r(l - v) 



<^^r = -T^ZTT^ sin(^((l - 2iy){lnr + Ko{Kr)) + - K2{Kr) }, (4.20) 



OziD = — 



27r(l - v) 

IJ.fl 
27r(l- 



cos^|(l-2i/)(lnr + i4:o(Kr)) + l-^+K2(Kr)|. (4.21) 



The fields (4.16)-(4.21) agree with the expressions given by Povstenko [20] in the framework of 

nonlocal elasticity if we use the identification k = 1/{tI). He used the two-dimensional nonlocal 
kernel (2.21) which is Green's function of the two-dimensional Helmholtz equation. Again, if 
one replaces flz by b, the stresses (4.16)-(4.21) agree with the stresses of an edge dislocation 
in cylindrical coordinates (compare with equations (3.27) (3.30) in [32]). The stresses (4.16)- 
(4.21) coincide with the classical ones far from the disclination core. The stresses (4.16)-(4.19) 
are zero at the disclination line. In principle we may discuss the extreme values of (4.16)- 
(4.21) in the xy-plane. For example, the stresses (4.16) and (4.18) have the extreme values: 
\arr\ ^ 0.260Kf§^^ at r ~ 1.494k-i and \a^^\ ~ 0.547K^§p^ at r ~ 0.996^-1. The 

stresses (4.20) and (4.21) have the values at r = 0: ct^^(O) ~ 2 "(i-iT) [(^ " 2iy){-f + ln f ) - i] and 

'^.^(0)^^n?^[(l-2t^)(7 + lnf)-i]. 
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(e) Kx 

Figure 1: The stress components of a twist disclination near the discUnation hne: (a) axx{0, y)i 
(b) a-yy{0, y), (c) (Jxy{x, 0) are given in units of hQ,zk/\^'k{1 — u)], (d) (Tzz{^, v) is given in units 
of fii}i^K/[TT{l — z^)] and (c) azy{x,0) is given in units of /if2/[27r(l — u)]. The dashed curves 
represent the classical stress components. 



The elastic strain is given in terms of stress functions 



, d'yj-uAf -dij -dyF 

27t 



Eij = —\ dlxf-y/^f dxF + d,g\. (4.22) 



-dyF dxF + d,g 
For the elastic strain of a twist disclination we find 

= -TTT^TTTT §1 (1 - 2-) + ^ + ^(.^ - 3a:^) (4.23) 



47r(l - v) 



~^{^-'') i^rKiinr) - ^(y' - 3a;2)i^2('*r)|, 



Eyy = il— ^l{i-2,y)-^-^(y^-3x^) (4.24) 

™ 47r(l-i/) r2\^ ' r2 ^2^4 V^^ J v ; 

- 2 - i.^ Atrial (Kr) + ^(y' - 3a;2)i^2(Kr)|, 

„ fl zx ( 2y'^ 4-99. .X 



47r(l — z^) r'^ \ K?r'^ 



2y2 2 

-V Kri^i(«:r) + — (x^ - 3y2) is:2 (Atr) 



Ezx 



= "MTn;) {<' - ^.-Xlnrt A-.C-)) + - ^4^1^ (2 - '''•■'-^^(.r))}^ (4,27) 

The components of the strain tensor have in the a:;y-plane the following extreme values (y = 
0.3): \Exx{0,v)\ ~ 0.308«4rffe;y at \y\ ~ 0.922k-i, |£;^^(0,y)| ~ O.OlOK^^il^ at |y| ~ 

0.218K-1, |£;j,j;(0,2/)| 0.054K4^i|^ at \y\ ~ 4.130k-i, and |£;,j,(x,0)| ~ 0.260«:4^;i^ at 
|a;| ~ 1.494k;~-^. It is interesting to note that Eyy{0,y) is much smaller than Exx{0,y) within 
the core region. The strain E^y has at r = the value: E^y (0) ~ 4,„(iL^) [(1 — 2i/) (7 -I- In |) — 5] . 
The strain (4.23)-(4.27) coincides with the result given by Gutkin and Aifantis [26-28]. The 
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dilatation Ekk reads 

In the xy-planc it has the extremum \Ekk{0, y) \ — 0-399^ ^^^^"^^'^^^^ at \y\ ~ 1.114k~^. The elastic 
strain can be rewritten in cylindrical coordinates as follows 



47r(l — v) r [ K'^r'^ J 



(4.29) 



E 



Q. z cos ip 



47r(l — u) r 



{i-;^ + 2^2m|, (4.30) 



= -TZTT—^ ( (1 - 2z^) + i - ^K^inr) - 2(1 - i^)nrKi{Kr) } , (4.31) 



47r(l - v) 

^ ~ 4n{l- 1^) sin<^|(l-2i^)(lnr + i^o(Kr-)) +^-i^2('«r)|, (4.32) 
^^^^~ 47r(l"-z.) cosv5|(l-2z/)(lnr + ifo(KO) +l-;^+^2(Kr)|. (4.33) 

The main feature of the solution given by (4.23) (4.33) is the absence of any singularities near 
the disclination line. This solution coincides with the classical ones far from the disclination 
core. 

Now we want to calculate the elastic bend-twist, torsion, disclination density and the rota- 
tion {ujz = —P[xy] s-iid LUy = P[xz]) of ^ twist disclination. They might be determined from the 
following conditions on the dislocation densities of the twist disclination: 

ay, = dxAf - dyw, = 0, (4.35) 

Q,, = i- ( + dl,g) - = 0, (4.36) 

- ^ {d^yF - dlyf) - = 0, (4.37) 

^ {dlyF + {dlJ - Z.A/)) - kyx = 0, (4.38) 

{dlyF + dl,g - a, {dlJ - t/A/)) - fc., = 0, (4.39) 



^vx — 



1 

2/1 



1 



Olzx = 



2^ {dl,F + dig + dlj) - kyy = 0, (4.40) 

-^d,{dxF + d,g)-kxz = 0, (4.41) 
1 

2/1 

= 2fc,i = 0. (4.43) 



azy = -w-d^zyF-ky, = 0, (4.42) 



The conditions (4.34)-(4.43) must be consistent with deWit's dislocation densities of a twist 
disclination. We will discuss this point in detail below. The Eqs. (4.34) and (4.35) look like the 
conditions for the dislocation density of an edge dislocation (see [32]) such that the elastic bend- 
twist kzx and kzy are compatible. Eqs. (4.41) and (4.42) are trivially satisfied. From (4.34)- 
(4.40) we may determine the elastic bend-twist. So we find for the non-vanishing components 
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of the elastic bend-twist tensor 



fcyx = - ^ ^ { 1 - AcrXi (kt)}, (4.44) 

fc.x = ^{(a;'-y')(l-«r-i^i('«r)) -KWifo(Kr-)}, (4.46) 

fc^y = ^ a:2/{2(l - KrXi(Kr)) - K\^Ko{Kr)}, (4.47) 

fc.. = - 1^ 4 1 1 - '^^^1 («0 } • (4-48) 



The shape of (4.44) and (4.45) is analogous to the elastic bend-twist of a wedge disclination 

given in [34]. It can be seen that (4.46) and (4.47) are singular at the disclination line r = 0. 
If one replaces flz by b, Eqs. (4.44) and (4.45) coincide with the elastic bend-twist of an edge 
dislocation (see [32,33]). The component (4.48) has no singularity at the disclination line. 

The clastic bend-twist tensor can be decomposed according to (2.8) into a gradient of the 
rotation vector and an incompatible part. We identify the incompatible part with the disclina- 
tion loop density. It is analogous to the decomposition of the elastic distortion of a dislocation 
into a gradient of the displacement vector and an incompatible distortion (see [30-32]). We 
find for the non-vanishing components of the rotation vector 

= ^ |<p(l - KrKi{Kr)) + ^ sign(y) Kr/^i («;r) | , (4.49) 

= - — ^ 1 1 - ktK^ {kt) I . (4.50) 

Here we use a single- valued discontinuous form for ip (see [13,22-24]). It is made unique by 
cutting the half-plane y = at x < and assuming ^p to jump from tt to — tt when crossing 
the cut. The far fields of the rotation vector (4.49) and (4.50) agree with deWit's expressions 
given in [13]. It yields signy = -1-1 for y > and aigny = —1 for y < 0. When y +0, the 
expression (4.49) is plotted in Fig. 2a. It can be seen that the Bessel function terms which appear 
in (4.49) lead to the symmetric smoothing of the rotation vector profile, in contrast to the abrupt 
jump occurring in the classical solution. It is interesting to note that the size of such a transition 
zone is approximately 12 /k which gives the value 6/k for the radius of the disclination core. The 
component (4.49) is discontinuous due to ip and (4.50) is continuous. The component (4.50) 
has in the xy-plane a maximum of Uz{x,0) ~ 0.399^2^2:/ [27r] at a; ~ — 1.114/k and a minimum 
of u)z{x,0) ~ — 0.399f2Kz/[27r] at x ~ 1.114/k and no singularity at the disclination core (see 
Fig. 2b). It can be seen that kzx, kzy and kzz are gradient terms of the rotation tJz- In 
performing the differentiations of the rotation Wy we obtain kyx and kyy plus excess terms which 
we identify as components of the disclination loop density. The non-vanishing components of 
the disclination loop density turn out to be 

'pIx = ^ n^xKQ{Kr) sign(2/)) , (4.51) 

fly = ^{/«'i/^s:o(Kr)(¥'- ^sign(j/)) +7r<5(y)(l-sign(a;)[l-KrJi:i(Kr)])}. (4.52) 

They contain the angle ip and the form is analogous to the plastic distortion of a disloca- 
tion (see [32]). Only the component ip*y has a ^-singularity at y = like the disclination loop 
density [13,35] ^ly = (17/2) (5(y)(l - sign(x)). 
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Finally, we find for the elastic distortion of the straight twist disclination 

/3.. = -rr^ ^((1 -^-) + K + 4-Ay'- ^^') (4-53) 

47r(l — v) r'^ [ K^r* 



47r(l — y) r'^ \^ r'' n'^r* 

- 2 (^(1 - I/) + KrJri(Kr) + ^(ar^ - 3y2)i^2M|, 

= -4^ - 2^) + ¥ + ^(^^ - ^^^) ^'-''^ 

/3 = JL^ £| / (1 _ 2z.) - - ^ (y2 _ 3^2N (4_5g^ 

47r(l - !/) r2 \^ ' r2 ^ ^ ' 

- ^ |</'(l " KrKiinr)) + ^ sign(y) KrKi^Kr)^ 

+ ^ |</'(l - KrKi{Kr)) + ^sign(t/) «;ri4'i(/tr)| , 
P.y = -^^^(^{(1 - 2.)(lnr + KoM) ^^^^(2 - «V^i^2M) |, (4.59) 

= -AMT^){^^ ~ 2^)(1^- + ^o(«r)) + ^ - ^^^(2 - «Vi^2(«r)) |. (4.60) 

Replacing Qz by 6, Eqs. (4.53)-(4.56) are analogous to the elastic distortion of an edge dislo- 
cation (see [32]). The components of the elastic distortion (4.57) and (4.58) contain the angle 

if in contrast to the dislocation case. But this is a typical property of a disclination. 
With Eq. (2.12) wc obtain for the effective Frank vector of the twist disclination 

fly{r) = j> [kyxAx + kyydy) = f2 |l — KrJCi (wr) | . (4.61) 

It differs appreciably from the constant value fl in the region from r = 0uptor~6/«; (see 
Fig. 3). In fact, we find ilj,(0) — and ^ly{oo) — f2. Thus, it is suggestive to take Tc — Q/k as the 
core radius of the disclination. The effective Frank vector Oy(r) of a straight twist disclination 
has the same form as the effective Prank vector Cl^ of a straight wedge disclination which is 
given in [34]. 

In the case of a twist disclination we obtain the following disclination torsion 

ax, = —zKo{Kr). (4.62) 

ZTT 

It looks like a dislocation density of a straight edge dislocation whose "Burgers vector" Q.z 
depends on the position z. At the point z = the dislocation density (4.62) is zero. The dislo- 
cation line of the edge dislocation coincides with the disclination line of the twist disclination. 
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-10 -5 5 10 

(b) KX 

Figure 2: Rotation vector of a twist disclination:(a) ojy{x,y +0)/f7, (b) ojz{x,0) is plotted in 
units of ilzK/[2TT]. The dashed curves represent the classical solution. 

Therefore, this dislocation density implies a dislocation line with changing Burgers vector in 
agreement with dcWit [13]. In the limit 1/k 0, dc Wit's classical expression a^z = ^zS{r) is 
restored. Consequently, we found that the straight twist disclination contains a certain amount 
of dislocation density (see also [13]). 

The elastic distortion gives rise to an effective Burgers vector 

(r) = {Pxxdx + p^ydy) = { 1 - KrKi (^r) } . (4.63) 

We see explicitly the changing of the Burgers vector on the dislocation and disclination line. 
The effective Burgers vector differs from the constant value J7z in the region from r = to 
r ~ 6/k. We find 6^(0) = and bx{oo) = Qz. In addition, the Burgers vector (4.63) depends 
on the position of z. At the position z = it is zero. From (4.61) and (4.63) we obtain the 
relation between the effective Burgers and Frank vector 

b^{r) = zflyir). (4.64) 

We find for the non- vanishing component of the disclination density (2.11) of a twist discli- 
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Figure 3: Effective Prank vector Q.y{r)/Q, (solid). 



nation 



(4.65) 



In the limit as k ^ ^ 0, the result (4.65) converts to the classical expression Qy^ = 5(r). 
The disclination density tensor (4.65) and the disclination torsion (4.62) of a twist disclination 
fulfil the Eq. (2.16) as follows 



(4.66) 



Since the dislocation density (4.62) and the disclination density (4.65) arc localized at the 
same position it seems that the dislocation density (disclination torsion) is coupled on the twist 
disclination. It is a characteristic quantity of a twist disclination which cannot be created by a 
pure edge dislocation without the presence of a twist disclination. In general, the disclination 
torsion is not independent of the disclination density of a straight twist disclination. Only in 
the xy-plane at ^ = the disclination torsion and the Burgers vector of the corresponding twist 
disclination are zero. 

If we use the decomposition (2.1) for the distortions (4.53)-(4.58), we may restore an effective 
displacement field and a properly incompatible distortion. The displacement is given by 

f TT 1 xy / 4 \ 1 

" 2^1^^^ ~ Kri^i(Kr)) + - sign(?/) KrKi{Kr) + _ ^ ' ^ + 2^2(Kr) J |, 

(4.67) 



^" 47r(l-i.) 

n 



|(l-2z.)(lnr + i^o(Kr)) + ^-^^^2^(2-«Vif2(Kr))|, (4.68) 



27r 



+ 



ip[l — KrKi{Kr)) + — sign(y) KrKi{Kr) 



{l-2iy){lnr-l + Ko{Kr)) - (2-KVK2(Kr)) 



2(1 - v) 



(4.69) 



These displacements (4.67)-(4.69) have no singularities at the disclination line. When y ^ 0, 
the Besscl function terms in (4.67) lead to the symmetric smoothing of the displacement profile 
in contrast to the abrupt jump occurring in the profile of the classical solution (see Fig. 4a). 
Eqs. (4.68) and (4.69) demonstrate the elimination of "classical" logarithmic singularities at the 
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Figure 4: Displacement vector of a twist disclination: (a) Ux{x, y +0)/[nz], (b) Uy{x, 0)/[ilz] 
with 1/ = 0.3 (c) Uz{x,y — > +0)/fl. The dashed curves represent the classical solution. 
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disclination line (see Fig. 4b where a constant term proportional to ln«; is neglected). When 
t/ — * 0, the Bessel function terms in (4.69) smooth the displacement profile in the core region 
(see Fig. 4c). It is interesting to note that the non-classical parts of the displacements (4.67) 
and (4.68) caused by the straight twist disclination with the Frank vector = (0, il, 0) coincide 
with the non-classical parts of the displacements due to the straight edge dislocation with the 
Burgers vector 6 = (6,0,0) if we replace Slz by b (compare with equations (3.47) and (3.52) 
in [32]). In addition, the displacement (4.69) coincides with the displacement —Uy of a wedge 
disclination (compare equation (43) for C = in [34]). The classical parts of (4.67)-(4.69) 
agree with the displacement given by deWit [13]. The displacement values should be detectable 
for nanoparticle containing twist disclinations. Therefore, one coiild compare the displace- 
ments (4.67)-(4.69) with experimental and simulated results. In performing the differentiations 
of the displacement (4.67)-(4.69) we obtain the total distortion f3fj = djUi. Using (2.1) and 
comparing the total distortion with the clastic one (4.53) (4.60) the excess terms of the total 
distortion may be identified with the plastic part. So the incompatible distortion can be found 
as 

= K^xKoinr) - | sign(?/)) , (4.70) 

= [K^yKoiKr) - | sign(2/)) +TrS{y)(l - sign(a;)[l - Kri^i (wr)] ) } , (4.71) 

= ^ K''x''Ko{Kr) - I sign(y)) , (4.72) 

P^y = --^ ^^KrKi{K.r) - K^xyKo{K,r)(^ip - ^ sign(y)^ - T:5{y)x{l - sign(a;)[l - nrKiinr)]^^. 

(4.73) 

Eqs. (4.70)-(4.73) satisfy the relation (2.7). The (5-terms in (4.71) and (4.73) have a similar form 
like deWit's plastic distortion [13] of a twist disclination = {Q.z/2) 5{y){l — sign(a;)) and 
j3^y = —{Q,x/2) 5{y){\ — sign(.T)). But now the singularity surface is not strictly bonded by the 
disclination line. The incompatible distortions (4.70) and (4.71) coincide with the incompatible 
distortion of an edge dislocation if we replace Q.z by b (see [32]) and (4.72) and (4.73) agree 
with the incompatible distortion —j3yx and —Pyy of a wedge disclination (see [34]). 
Using (2.15) we obtain for (4.70)-(4.73) the following decomposition 





Pxx — 


ZWyx, 








Pxy — 


ZWyy, 








Pzx = 


-XWyx, 








hy = 


(i>Zy - XWyy, 




(4.74) 


into the translational gauge field 












(t>zy = - 


~ KrKi{Kr), 

ZTT 




(4.75) 


and the rotational gauge field 










Wyx 


= -Vyx 


Wyy = - 


'V*yy 


(4.76) 



Thus, the negative disclination loop density (4.51) is equivalent to the rotational gauge poten- 
tial (4.76). 

5 Conclusion 

The field theory of elastoplasticity has been employed on the consideration of a straight twist 
disclination in an infinitely extended body. We were able to calculate the elastic and plastic 
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fields. We found that the elastic stress, elastic strain, elastic bend-twist, dislocation density 
and disclination density are continuous and the displacement, plastic distortion, rotation and 
the disclination loop density of the twist disclination are discontinuous fields. Exact analytical 
solutions for all characteristic field quantities of a twist disclination have been reported which 
demonstrate the elimination of "classical" singularities at the disclination line. The disclination 
core appears naturally as a result of the smoothing of the rotation vector profile. In addition, 
wo pointed out and discussed the relation between the twist disclination with Frank vector 
CI = (0, f2, 0) and an edge dislocation with Burgers vector b = (6, 0, 0). We were able to calculate 
the effective Frank and Burgers vector of the twist disclination. The force stress of a twist 
disclination calculated in the field theory of elastoplasticity agrees with the stress calculated 
within the theory of nonlocal elasticity and strain gradient elasticity. The reason is that in all 
three theories the fundamental equation for the force stress has the form of an inhomogeneous 
Helmholtz equation (see Eq. (2.19)). The solutions of a twist disclination considered in this 
paper could be help in studies of mechanical behaviour of nano-objects including nanotubes 
and nanomembranes and of disclinated nanoparticles. Last but not least, using the geometric 
framework of /5'0(3)-gauge theory of defects we have found the translational and rotational 
gauge fields of a twist disclination. It turned out that the (negative) disclination loop density is 
equivalent to the rotational gauge field. In general, the (negative) gauge fields may be considered 
as the plastic parts in the field theory of elastoplasticity. 
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